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A field theory to describe electrostatic interactions in soap films, described by electric multi-layers 
with a generalized thermodynamic surface-charging mechanism, is studied. In the limit where the 
electrostatic interactions are weak this theory is exactly soluble. The theory incorporates in a 
consistent way, the surface-charging mechanism and the fiuctuations in the electrostatic field which 
correspond to the zero frequency component of the van der Waals force. It is shown that these terms 
lead to a Casimir like attraction which can be sufficiently large to explain the transition between 
the Common Black Film to a Newton Black Film. 



I. INTRODUCTION 

Soap films are naturally occurring examples of diffuse double layers |^ . They consist of two surfaces formed by 
hydrophobic surfactants which accumulate at the air water interface due to the hydrophobic nature of the hydrocarbon 
chains of the surfactant molecule. In the classic experimental set up to measure the disjoining pressure in soap films, 
the film is placed in a cell in contact with a bulk solution, usually in a porous frit ^, Q]. A capillary connecting 
the bulk solution to the outside of the cell allows the pressure in the cell to be varied and a direct measurement of 
the disjoining pressure Pd, generated by interactions within the film, is thus possible. The thickness of the film as a 
function of Pd, along with certain aspects of its structure such as the electronic density, can then be measured via 
X-ray or optical methods. The interface formed by the surfactant can lead to the formation of a surface charge and, 
in addition to the surfactant, the bulk can contain an electrolyte such as salt. 

The simplest model of a soap film is of two charged surfaces separated by an electrolytic solution. Electrostatic 
interactions therefore play an fundamental role and can be taken into account at the simplest level of approximation 
by the mean field Poisson Boltzmann equation ||. In the simplest models, the surface charge or surface potential 
is taken to be fixed and is hence a fitting parameter of the theory. More elaborate theories take into account the 
mechanism of surface charging 1^, ^, 1). Steric effects between the counter-ions in the vicinity of the surface, known as 
the Stern layer, may also be incorporated |^. In general the disjoining pressure predicted by such theories is repulsive. 
In a detailed experimental study of soap films formed by ionic surfactants such as sodium dodecyl sulphate (SDS) 
[^, it was shown that a modified Poisson Boltzmann theory incorporating an energetic surface-charging mechanism, 
via an attractive free energy potential at the film surface, can predict the surface tension of the bulk solution. Then, 
with no free fitting parameters remaining, the theory predicts the disjoining pressure of the films made from this 
bulk, just up to a pressure where the Common Black Film (CBF) collapses and forms a Newton Black Film (NBF). 
The NBF is an extremely thin film where experiments show that there is an extremely small separation between 
the two surfactant surfaces (about one layer of water molecules across). The transition between the CBF-NBF 
is experimentally interpreted as a first order phase transition [Q, and the body of theoretical work supports this 
interpretation. The rough picture, coming from the approximate theories that exist, is that at large distances the 
disjoining pressure is repulsive and stabilizes the CBF. At closer inter-surface distances attractive van der Waals forces 
some into play which are responsible for the eventual collapse to the NBF. 

The aim of this paper is to analyze the role of electrostatic fluctuations in the presence of generic surface-charging 
mechanisms and spatially varying dielectric constants in diffuse layer systems. Previous studies incorporating the 
fluctuations about the mean field Poisson Boltzmann theory in simpler models have revealed that these fluctuations 
lead to attractive interactions |l2[ The analysis here is adapted to the limit of weak electrostatic interactions. 
We use a field theoretic formulation of the problem, which in addition to producing new results for fluctuating surface 
charges, allows one to recover, in a powerful and universal way, many results already established in the literature. In 
addition our formulation is well adapted to develop a perturbative expansion which allows one to incorporate higher 
order interactions. 

The Poisson Boltzmann theory may be supplemented by including the contribution of dispersion forces or van der 
Waals forces calculated from the Hamaker theory based on pairwisc dipole interactions and, at a more sophisticated 
level, via the continuum Lifshitz theory . The resulting DLVO [Q, ^ theory has been very successful in describing 
the physics of electric double layers. However, the splitting of the overall interaction into a Poisson Boltzmann static 



interaction and the van der Waals forces as two distinct and independent interactions is artificial from a global view 
point. The fluctuations of the full theory of the electrostatic interactions generate the zero frequency van der Waals 
interaction. However since the ion distributions do not respond to high frequency fields, the non zero frequency 
contributions to the van der Waals forces can be taken to be effectively independent of the ionic distribution. The 
zero frequency contributions between neutral surfaces or surfaces with absorbed mobile ions or dipoles can be treated 
and one finds that the zero frequency van der Waals forces become screened by the presence of an electrolyte . 
In the context of the Debye-Hiickel approximation these van der Waals forces have a formally identical origin to the 
Casimir effect, where the suppression of fluctuations of the electromagnetic field due to the presence of two surfaces 
leads to a net attractive force. The appearance of such generalized Casimir forces in soft condensed matter systems 
and many other contexts is now well established ItI p^. 

The theory proposed for electrolytic soap films in 1 5[ can in fact be solved exactly in one dimension , via path 
integral techniques originally developed in pc| ]. In this exact theory one finds that the disjoining pressure is repulsive 
at large distances but attractive contributions come into play at smaller inter-surface separations and can lead to a 
disjoining pressure isotherm predicting a collapse to a thinner film state at a certain thickness, as is seen in experiments. 
Attractive forces in the two dimensional form of this model have also been found l2l|| at a particular critical temperature 
where the model is exactly soluble. In this paper we revisit the theory of ||^ but in three dimensions. We use a field 
theory representation of the system in the weak coupling limit, equivalent to the region where the Debye-Hiickel 
approximation is valid, where the field theory is free and one may decouple the Fourier components of the field and 
apply standard path integral results. This allows one to incorporate a thermodynamic or energetic surface-charging 
mechanism in a straightforward way and also allows one to take into account spatially varying dielectric constants. 
In principle the solution of the free field theory can be written down in terms of a functional determinant which may 
be evaluated via functional techniques |2^, |2^, |2^; these functional techniques could also be applied here, 

however the path integral method gives a very compact and rapid solution to the problem. 

Recent experimental studies of the CBF-NBF transition have been carried out on nonionic soap films made with 
the surfactant C12EQ Q. These films can be stabilized in the presence of small concentrations of electrolyte and are 
thus in a region where the linearized Poisson-Boltzmann equation can be used for the mean field treatment and also 
where the free or Gaussian field theory we examine should be valid. 



II. THE MODEL AND SURFACE-CHARGING MECHANISMS 

Here we explain the derivation of the type of model proposed in ||^, |l^ but with a generic surface-charging mecha- 
nism. We consider the case of a monovalent electrolyte for simplicity. In the presence of the surface, a surface charge 
may be generated by two basic mechanisms. Firstly, as was taken to be the case in there may be an affinity 

for one of the species to be on the surface. In the case of SDS, or other ionic films, this species is the soap tail as, 
due to the hydrophobic effect, the hydrocarbon tail lowers the free energy of the system by leaving the aqueous core 
and entering the air environment at the exterior of the film. The head of the surfactant, SO^ for SDS, is however 
negatively charged, thus leading to a surface charge. In the case of nonionic soap films the surface-charging mechanism 
is due to a difference in mobility of the ions. Hydration effects can make one species effectively larger than another 
thus leading to a steric repulsion at the surface and hence an effective charge is induced by the unneutralized presence 
of the other species. The enhancement of the repulsion between the surfaces of ionic Aerosol-OT films due to more 
hydrated counter-ions, e.g. LiCl instead of CsCl, has been experimentally demonstrated in ||2^. Let us consider a 
model with one salt species (water is also partially dissociated and indeed plays an important role in the absence of 
salt), the generalization to several species is straightforward. 

Let us denote by Vi{x) and V-i{x) the effective potential for the cations and anions respectively due to the surfaces 
of the film. The Hamiltonian for the electrolyte system is thus 



H=-eY^^{^M^ + Y.^,A^^) (1) 

i i 

where the first term is the electrostatic energy with potential t/'j — il if particle i is a cation/anion respectively 
and Xi its position. The second term is the boundary interaction. We denote by L the perpendicular distance between 
the two surfaces (in the direction z) and by A the surface area of the film (in the plane (x, y)). We denote hy Ax L 
the region inside the film and by A x T the region outside the film. The total length of the system in the z direction 
is denoted by U and therefore T — U — L. To start with, we ignore the fact that just outside the surface one has a 
density of hydrocarbon tails (this will be taken into account very simply later on). The calculation for a triple layer 
in the absence of surface charging and electrolyte can be found in |2^, In this two layer picture therefore, the 
dielectric constant has the form e(x) = e for x G ^ x L where e is the dielectric constant of the electrolyte solution. 
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which in the dilute hmit we shall take to be the dielectric concentration of water. Outside the film we have e(x) = cq 
(x S A X T), where eo is the dielectric constant of air. The difference between the internal and external dielectric 
constants means that one needs to take into account image charges and the zero frequency terms of the Lifshitz theory, 
however this is automatically incorporated in the field theoretic formulation of the theory. 
The electrostatic potential ip satisfies the Poisson equation 



V-e(x)V7/'(x) ^-eyg,5(x-x0 (2) 



We follow the standard method for converting a Coulomb system to a Sine-Gordon like field theory by performing 
a Hubbard-Stratonovich transformation with an auxiliary field (p to obtain the grand canonical partition function 



d[<i>]e^piSm (3) 



where 



S[(j)]=-^ j /3e(x)(V0)2dx + fi [ exp{i(3e(j)- (3Vi{x))dx 

^ J{T+L)xA J LxA 

Ai / exp(-ie/30-/3y_i(x))dx (4) 

J LxA 

and where is the fugacity of the cations and anions. It should be noticed that the Hubbard Stratonovich- 
transformation is over all space and hence the first, kinetic, term of the action is an integral over all space denoted 
(T + L) X A and the second two terms are the interaction terms restricted to the film region LxA. We notice that 
the functional integral gives the ion/ion interaction upon performing the integral, but in addition there is a term 

j d[(j)\ exp j ■ e(x)V0dx^ = (det(-V • e(x)V))~* (5) 

This term is however the zero frequency contribution coming from the Lifshitz theory and in fact should be there. 
This term naturally arises when one considers the full Quantum Electrodynamics (QED) of the system. If one fixes 
the positions of the ions and ignores the magnetic part of the Lagrangian, the time independent (zero Matsubara 
frequency) purely electrostatic part of the Lagrangian is 

C[iP] = e(x)(VV')2dx-e^(?,^(x,;) (6) 

hence the thermal field theory for the field has a partition function HJ] 

Z = I d[V']exp(/3£[V']) (7) 



If one now takes a classical trace over the ion positions and passes to the grand canonical ensemble, one finds the 
expression for the grand potential above after changing the axis of the functional integration via ^ — *■ —icj). If one 
wishes to take into account the non zero frequency Lifshitz terms one proceeds as above but keeping the full QED 
action and introducing the bosonic Matsubura frequencies ujn — 2nn/f3?i ||3l|. Hence the grand partition function S 
contains the ionic interactions and zero frequency van der Waals contributions in the system. We define a Stern layer 
to be the region where either or both of the potentials V±i are non zero and take the width of this region to be 6. 
The action then has the form 



S[(t)]^ - \ j (3e{x){V(t)fdx + 2n [ cos (/3e0) dx 

2 J{T+L)xA JlxA 



ji j exp (i/3e0) (exp (— /3Vi(x)) — 1) dx + /i / exp [ifiecf) (exp (— /3Vi(x)) — 1) dx 

J[0,i5]xA J[L-S,L]xA 

/ exp(-i^e(/)) (exp(-/3y_i(x)) - l)dx + ^ / exp (-i/5e^) (exp (-^y_i(x)) - 1) dx (8) 

J[Q,S]xA J[L-5,L]xA 
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We now take 6 to be small and take the limit 5 — > choosing 

(exp (-/3y±i(x)) - 1) ^ i^UHz) + S{L - z)) (9) 

where S{z) is the one dimensional Dirac delta function. Clearly /j.^. is positive/negative for V±i negative (attraction 
of the species towards the surface)/ positive (repulsion of the species from the surface). Note that dimensionally 
[a*±] = [m] [^]- III this simplified limit, the action is now 



S[(l)]=-^ [ P€{x){V(l)fdx + 2fi I cos {(3e(l))dx 

2 J{T+L)xA JlxA 

+ {5{z) + 6{L - z)) exp (i/3e0) dx 

JlxA 

+ H*_ I ((5(z) + (5(L - 2;))exp(-«/3e^) dx 

JLxA 



(10) 



In the weak coupling limit we proceed by expanding the action to quadratic order in all terms, yielding the Gaussian 
action 



Sg [<t>] = -l f pe{z){V<j>fdx-l f pem^cp^dx 

+ iX it (j) dxdy + (j) dxdy] - ]-(3em'^'y ( [ (j)^ dxdy + [ (j)^ dxdy] + 2fiAL + 2{fi*_^ + fi*_)A (11) 

\Jz=0 Jz=L ) 2 \Jz=0 Jz=L ) 



Here m = ^J2iie'^j3/e corresponds to the Debye mass (the inverse Debye length) in the dilute, weak coupling limit, as 
in this limit /U ~ p, where p is the bulk electrolyte concentration. In addition we have 



and we recall e(z) = e for e [0, L] and is equal to eo elsewhere. We now express (j) in terms of its Fourier decomposition 
in the r = {x, y) plane. 

<^(r,^) = -^^^(p,^) (13) 

where \i A = ly.l imposing periodic boundary conditions yields p G ^{fixjUy) with {nx,ny) & Z^. In terms of the 
fields (p, the action decouples and one obtains 

Sg = So[4>m + E ^p[^(P)] + + 2(m; + (14) 



where 



and 



Soim = -I i^^^^^j dz-lpem^ Jj\0,z)dz 

+ iXVA (0(0, 0) + ^(0, L)) - i/3e7m' (j>''{0, 0) + ^^(o, L)^ (15) 

^■'[^(P)^ = -I Lj<'^ +p-^(p,.)^(-p,.)) dz 

- If^em^-f [ 4>{p,z)^{-p,z) dz (16) 
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The system therefore decomposes into a system of simple harmonic oscillators and 

S = AAexp {2iiAL + 2A(/x; + ^^1)) f d[0(O)] expi^o] '^[^(p)] exp[5p] (17) 

where in the following J\f will be used to note L independent normalization factors. Each oscillator (labeled by p) 
has a time dependent Hamiltonian 

i/p = -^^ + >/(.W.)^^ (18) 



where M{z) = Pe{z) and w(p, z) = y^p^ + rn^ for z E [0, L] and uj{p, z) = |p| = p for z ^ [0, L]. Here X represents 
the field 0(p) and z corresponds to the temporal coordinate. The mode p = is slightly more complicated due to 
the presence of a linear source term. One may write 

d[0(O)]exp[5o] = 

A^oTr exp (-([/ - L)He) exp [iX\fAX - \l3em^^xA exp {-LHp) exp {i\\fAX - ^/3em27X2 



(19) 

where Hp is the Hamiltonian outside the film with mass Me = Peo and loe = and Hp is the Hamiltonian in the 
film with Mp = (3e and top = m (subscripts E and F will be used to denote the masses and frequencies external and 
internal to the film respectively) . In the limit where U oo only the eigenstate of Hp of lowest energy survives the 
thermodynamic limit (we shall see later that this is compatible with the constraint of electro-neutrality within the 
film at mean field level). We recall that the ground-state wave function of the simple harmonic oscillator is given up 
to a normalization by 

, MajX"^^ , , 

MX,M,u;)^eM —) (20) 

and the corresponding ground-state energy is ■^ui. In the case where p = one has the ground-state wave function 
'iljoiX,ME,0) = 1. Using this in Eq. (|l|) thus yields 



d[0(O)]exp[5o] = 

TVo J exp (^iXVAX - ^Mpm^-fX^^ K{X, Y, L, ujp, Mp) exp (^X^/AY - ]^Mm'^-iY^^ dXdY (21) 
where K{X, Y.L,uj, M) is the Feynman kernel for the simple Harmonic oscillator pO[ | and is given by 

K(X,Y,L,uj,M)^ ( --MLJC0th(i^L)(X^ + Y^ -2XYsech(ujL))\ (22) 

\27rsmh(cjL) J \ 2 J 

Performing the above integral yields a term extensive in A coming from the linear source term and a non extensive 
fiuctuation term, which we shall absorb into the normalization term, yielding 



In 



(/d[^(0)]exp[5o: 



A^(m7sinh(mL) -|- cosh(r7iL) -f 1) 

(23) 



f3em{{m?^'^ + 1) sinh(mL) + 2m7C0sh(mL)) 



This term can be simplified slightly giving the zero momentum contribution to the grand potential per area of film. 
Writing J = - ln(S)/(/3^) = Jc + Ep ^, where 

J, = -|(Mi + /i; + M;) (24) 



is ideal part of the grand potential, we find 

sinh (^) + m7Cosh {^) 



= , (25) 
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In fact we shall see later, as should be expected, that this term is the mean field contribution to the theory. 

The non zero momentum modes only contribute to the fluctuations but summed together give an extensive contri- 
bution. The contribution to the grand partition function of one of these modes is 



d[(^(p)]exp[S'p 
KTr 



exp(-— ^^^1^^) exp {~{U - L)He) exp 



exp (-LH^^) exp ( -i/Jem^X^) exp(-^^^^^^) 



(26) 



Note that there is no under-counting in this contribution as the field cj) was real and hence ipiji) — (f>{—p). One thus 
obtains 



d[(t>{p)] exp[5'p] = A^p exp 



U-L 



-UJe X 



exp 



Meu^eX' l_Mp^2^xA K{X,Y,L,LUE,ME)e^v ( MEi^EY^ ^-Mpm^^Y^ ) dXdY (27) 



2 2^ 
This yields an L dependent contribution to the grand potential per unit area of 



exp 



(28) 



The total contribution from the modes p 7^ is the grand potential associated with the static electrostatic fluctuations, 
or the zero frequency van der Waals force, hence we write J^dw = Sp^^to "^p- should remark however, given the 
unifield treatment used here, this notation is slightly arbitary as with the definition used here Jvdw contains terms 
coming from the ionic fluctuations. The first term in the expression ( p8|) will clearly lead to an ultraviolet divergence 
in the grand potential, however the disjoining pressure Pd{L) is given by the difference between the film and bulk 
pressures 



Ik 



dJ{L) 



8L 'ii-^(^)/^ dL 
After summing over all values of p, by passing to the continuum for large we obtain the divergence free result 



dJ*{L) 



(29) 



eaP 



+ em?"i — t\J p^ + rn^ 



eqP + em'^'-f + ey^p^ + m^ 



cxp(— 2L-\/p2 



2771/i7 

sinh I 



cosh (- 



^) + m7Cosh I 



(30) 



From Eq. ( |29| ) it is clear that J*{L) is the effective interaction potential between the two surfaces. 

The bulk grand potential contains the extensive contribution which itself contains the divergence. By taking the 
difference we have concentrated on the part relevant to the calculation of the disjoining pressure. However, it is 
instructive to remark that the extensive term in Eq. ( p8| ) is the term in the bracket that is explicitly proportional 
to L. It can be shown that on integration this term does indeed give the standard Debye-Hiickel expression for the 
pressure. It is, in fact, the one- loop contribution similar to that explicitly discussed in the ID model in The 
difference in the 3D case is that the integral is divergent and so must be renormalized. The divergence for this integral 
is due to the self-energy of the charge distribution included in the original definition of the partition function. It can 
be removed by imposing the renormalization condition the the fugacity /i be chosen so that the density is given by 



P = - ^ 



11 d log ^ 



/i is then a bare quantity which is a function of the Ultra- Violet cutoff. It can be shown that when the appropriate 
one-loop calculations are performed the bulk pressure is finite and is given by 



/3Pbuik 



P 



m 
24^ 
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Where they are important the UV divergences will ultimately be regulated by the finite size of the ions. However, 
these contributions will be corrections and it will be possible to estimate their importance using perturbative methods. 
The details will be postponed to a further paper which will report on the perturbative approach in general. The 
renormalization of /i is not of relevance to what follows but is central to higher-order corrections. 

III. VARIOUS LIMITS 

To demonstrate the generality of the formalism here and its compatibility with well established results in the 
literature derived by other means, we shall consider various limits of the formulas derived above. 

A. Mean Field Limit 

The mean field theory for the action is obtained from the equation 

S 



SG[<f>] = (31) 

This yields the equation 

/3V • eV(/) - m^Pecl) + iX{S{z) + d{z - L)) - I3em^-i(l){5{z) + 5{z - L)) = (32) 
within the film. Outside the film one has 

/3V • eoV0 = (33) 

As is usual for Euclidean Sine-Gordon field theories, the physical saddle point is imaginary, one writes <f> = iijj and 

it turns out that V' is the mean field electrostatic potential, this is also clear from the static QED formulation (|^). 



In addition the mean field solution depends only on the coordinate z. From Eq. (33) one finds that outside the film 
dip / dz = which is simply the condition of electro-neutrality in the film. Inside the film (strictly away from z = 
and z — L) one has therefore 

fie (^^V - = (34) 

which is simply the linearized Poisson Boltzmann equation. The solution to this equation which gives a potential 
symmetric about the films midpoint L/2 is 

il> = Ccosh ^rn(a; — (35) 

Integrating the mean field equation between 2 = 0~ and 2 = 0^, and using the condition of electro- neutrality, one 
finds 

- U=o+ - Pem'^iPiO) + A = (36) 



This allows one to solve for C giving 



Substitutin g th is mean field solution into the expression (^T|) for Sq, yields the mean field action which is exactly the 
expression ( P5| ) obtained from the zero momentum contribution in the previous analysis. 



/3em(sinh(2|^^) + m7COsh(^)) ^ ' 



B. No electrolyte no surface-charging 

This amounts to calculating the zero van der Waals forces across a slab of dielectric constant e while the dielectric 
constant of the external media is given by eg. Here one has 
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7 rj~i rOC f / \ 2 



J*^^) = ^ 7^ dp pin I 1 - (^-^ ) exp(-2Lp) ) (38) 



The integral above is easily evaluated by expanding the logarithm leading to 



n— 1 ^ 

This yields the attractive, Casimir like, disjoining pressure |lq, 05 



2n 



n— 1 ^ ^ 

C. No surface-charging 

If one considers the limit where there is no interaction between the ions and the surface one has = = and 
hence 7 = 0. This is the case of an electrolyte confined by two neutral surfaces (in both the electrostatic and chemical 
sense) but with a medium of dielectric constant eo out side the film and dielectric constant e inside the film. One thus 
obtains 



rW^'-f rdpplJl- h^-'^^l^ ] exp(-2LV^^) (41) 

This thus recovers the result of [ p^ demonstrating the screening of the zero frequency van der Waals interaction in 
the presence of an electrolyte. 

D. Localized Surface Charge Fluctuations 



It is easy to see that the term proportional to 7 in the action ( |11[) is a term proportional to the surface charge 
fluctuation about its mean value. The effect of surface charge fluctuations in the limit where the counterions are 
highly localized about the surface was considered in |^ and In |2|] the problem was considered in the strong 
and weak coupling Debye-Hiickel regime, in this regime the electrostatic fluctuations lead to the computation of a 
mathematically identical determinant in both |2q] and [p and the electrostatic fluctuations in our model also have the 



same mathematical structure. The differences between |6, 25 and our model are in terms of the physical interpretation 



at this level. In the case of non varying dielectric constant (for instance the case of two membranes interacting in 
water) the contribution to the disjoining pressure coming from the electrostatic fluctuations or the zero frequency van 
der Waals term p^j""^^^ (2^) is given by taking m = within the film. Immediately one obtains 

27r J [exp(2pL)(l +pA)'^ - IJ 



where A = 2/m^7. The expression ( [4^ ) is up to the definition of A the same as obtained in and ||2^ and g]. 

In the problem of localized surface-charging but with varying dielectric constant, that outside the film being 
different to that within the film, was also considered. In |2^ the model of |^ was studied incorporating the same 
dielectric constant variation. However at a mathematical level the calculation of the electrostatic fluctuations in 
|^5|,^^ and here is again identical. Here we find that the contribution from the electrostatic fluctuations is 

^r"^ w = / dp. — (43) 



2tt 



exp(2px) j^::::^ ) -i 



This equation agrees with the result of (in the Debyc Hiicklc regime) which is obtained there by calculating the 
contribution of the ionic fluctuations and then adding the contribution coming from the zero frequency Lifshitz terms. 
In psf only the ionic fluctuations were included and hence the result differs from ours and that of psf by the zero 
frequency Lifshitz term. 
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E. Including the Non-zero Frequency van der Waals Contribution 

Let us mention that if one wanted to now take into account non-zero frequency contributions from the Lifshitz 
theory, which are in principle uncoupled or very weakly coupled to the static ionic fluctuations one should use a total 
potential 



Jt{L) = J*{L) + Ju{{w, L) (44) 

where JLif (w, L) is the contribution to the dispersion potential coming from the Matsubara frequency w in the absence 
of electrolyte. As mentioned previously, in the formulation here, the zero frequency Lifshitz and ionic components are 
treated together in the grand partition function S. 



IV. THE TRIPLE LAYER 



As mentioned earlier, the change in the dielectric constant due the the presence of the hydrocarbon layer formed 
just outside the surface of the soap film by the surfactant should also be taken into account. If one considers this to be 
a region of size h at either interface, but inaccessible to the aqueous solution one now has an different Hamiltonian Hh 
(where the subscript H denotes the hydrocarbon region) in these regions characterized, for a mode p, by uir = P and 
Mh = /3ei, where ei is the dielectric constant of this hydrocarbon region (strictly this is a mixture of air/hydrocarbon 
chains). Using the formalism developed here one finds that the mode p yields a contribution 

d[4>{p)]cxp[Sp] = A/pexp (^-{- — ^ — —)ujE^ X 



2 K{X, Y, h, LOH, Mh) exp i--Mpm'^Y- 



/exp(- 

K{Y, Z, L, LOF, Mf) exp ^-^Mfto^tZ^^ K{Z, W, h, uh, Mh) exp (^_ ^eujeW ^ ^xdYdZdW (45) 

to the grand partition function. 

The p = mode or mean field contribution is unchanged by the addition of the hydrocarbon layer and after some 
algebra on finds that 



where 



with 



mL) = ^ fc^ppln l- ^^(^)^ + <-^^-vglg) ) exp(-2Lv^^) 

cosh (2^) , ^ 

+ '^muf T-r^ — i—TT (46) 

sinh(2|^^) +m7cosh(2^) ^ ' 



_ l + Aexp(-2p/i) 
^(^^ - 1-Aexp(-2M) (^^^ 



A = f^ (48) 
eo + ei ^ ' 

Note that if eo = ei we recover the double layer result obtained earlier. In the limit /i ^ we obtain eiB{p) eo 
and when /i — > oo one has eiB{p) a as it should. A key, physically illuminating step, in the algebra mentioned 
above is the evaluation of the term 

nY) = j eM-\MEOJEX'')K{X,Y,h,MH,oJH)dX (49) 
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One finds that 



exT) (- iy2 ( Me^eMhujh coth(cjgfe) + M^ul \ \ 

V 2 V MEUJE + MHUJHCoih{ujHh) )) ^ ' 

The normaHzation here is unimportant for the calculation of the disjoining pressure (as h is taken to be fixed). We 
can now interpret '0* as the new effective ground-state wave function associated to the mode p entering the film. As 
there is no electrolyte in the hydrocarbon region one has that loh = uje ~ P (the case where electrolyte was present 
could be easily handled by the same formalism) and hence one finds that the effective or renormalized ground-state 
wave function entering the film is 

Clearly this leads to a renormalization of Mh with respect to what it would be if the region of hydrocarbons was in 
fact of infinite size, and consequently one finds a p dependent renormalization of the dielectric constant ei as 

eo coth(pft,) + ei \ 
eo + ei cot 



e,(rt = e,B(rt = .i(-^^^^j (52) 



which is easily seen to be in agreement with the definition of B{p) above. This calculation can easily be generalized 
to an arbitrary number of electrolyte free layers simply by calculating recursively the value of B(p) generated by 
the composition of all the layers before arriving at the penultimate layer before the film. In the case of no surface 
interaction and no electrolyte this result agrees with the calculation of the zero frequency contribution to the van der 
Waals force in |, ||. 

In what concerns the the momentum dependence of e^{p) one finds that 

e^{p) w eo for i > /i (53) 
P 

e^{p) « ei for i < /i (54) 
P 

Hence the long wave length modes entering the film behave as if the hydrocarbon layer was not there and the short 
wave length modes behave as if there was only the hydrocarbon layer present. 

V. BEHAVIOR OF THE DISJOINING PRESSURE 

In this section we shall examine how the disjoining pressure is efi^ected by changing the various physical parameters 
of the theory. We write 

p.=pr+pr^) (55) 

where the first term is the mean field or zero mode contribution P^^^ — —dJo/dL and pj"'''^' — —^J*^^,/^L is the 
zero frequency van der Waals contribution coming from the modes p such that |p| ^ 0. 

A. The Mean Field Contribution 

The mean field contribution to the disjoining pressure is 

^ ' /3(sinh(2|^i)-Hm7COsh(^))2 ^ ' 
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which is clearly always positive. This repulsive component is generated by an effective surface charge which is present 
when fi*^ — 7^ 0, that is to say there is an asymmetry between the cation/anion afhnities or repulsions at the 
surface. At large inter-surface separations one finds 

which has the standard Poisson-Boltzmann exponential decay with the characteristic length scale of the Debye length 
Id = \/m. If 7 is negative 7 = —7' (with 7' > 0) the magnitude of the repulsion is enhanced. Here there is a critical 
value Lc where the mean field component to the disjoining pressure diverges. One finds that 

L.^UJ\±^\ (58) 



TO V 1 — "I7' 

The divergence here is not physical as it can be avoided by keeping higher order terms in the expansion of the surface 
terms in the full action ( p^ ) . However it does indicate an enhancement of the repulsion due to the presence of a Stern 
layer of effective length 5' — 7' from the considerations in section (II.). If the theory were applicable for large L and 
Lc were small, then one finds 

Lc w 25' (59) 

that is to say that the pressure should rise rapidly when the two effective Stern layers come into contact, thus giving 
a large repulsive term in the disjoining pressure before the surfactant surfaces actually come into contact. This image 
is valid in the range where mS' ^ 1, that is to say the width of the Stern layer is much smaller that the Debye length. 
It is clear that the limit taken in Eq. (^ is only valid in this case. In the case where m5' > 1 one must treat the 
Stern layer as a continuum and introduce two new layers of finite thickness. 

In the case 7 > the magnitude of the repulsive part of the disjoining pressure is decreased as 7 increases. At 
L = one finds the disjoining pressure 

Pf (0) = MbT(^^)' (60) 

Therefore at small interface separations one needs to incorporate steric repulsion between the surfactants in the 
two surfaces to prevent a collapse to a zero film thickness. 

B. The Zero Frequency Van der Waals Contribution 

As mentioned previously, the presence of an electrolyte in the film leads to a screening of the zero frequency van der 
Waals interaction and simple expressions for pj"''^) as a function of L do not exist due to the presence of a second 
length scale the Debye length (= 1/m). For simplicity we consider just the bilayer model. In the case L 3> 1/to, the 
disjoining pressure is dominated by modes such that p ^ 1 and one obtains to leading order 

jivdw) ^ _ ksTm^ / TO7 - 1 ^ ^ 
4:TTmL \ m7 + 1 



p,uay.> .^^^^ . exp(-2TOL). (61) 



We see that the prefactor controlling the strength of the long-distance attraction depends crucially on the value of 7 
and hence the surface-charging mechanism. The exponential decay is however twice as rapid as that of the mean field 
contribution, meaning that in the thick film regime the zero frequency van der Waals attraction is dominated by the 
mean field repulsion term. 

In the limit L ^ 1/to one finds 

n—l ^ ' 

This expression is independent of 7 and consequently coincides with that given by in jl5j . Hence we see that it is only 
in the regime of very thin films that the zero frequency van der Waals force takes the Hamaker form pj"'"^' 
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VI. ORDERS OF MAGNITUDE AND COMPARISON WITH EXPERIMENTAL DATA 



The theory presented so far is the hnearized version which is equivalent to the Gaussian approximation and corre- 
sponds to the free field formulation. We should expect that a detailed comparison with experiment will not be wholly 
successful since we have yet to include non-linear effects. However, it is instructive to show qualitatively how well the 
model performs and to this end we study one case for which the parameters are typical of experiments. We consider 
a soap film in air consisting of ionic liquid bounded by two thin hydrocarbon layers of hydrophobic surfactant. We 
idealize the model to be a bilayer system where we do not account for the non-zero thickness of the hydrocarbon 
layers. The liquid is composed of water with a dissolved salt, such as NaCl. The hydrophobic surfactant is of one 
charge only and so we set /i!^ = fi* , pL*_ = 0. The fugacity /i of the cations and anions is chosen by the experimenter 
and determines the Debye mass m through m? — 2iie^ /kT. For the free energy we use Eq. ( p8|) which we recast in 
terms of dimensionless variables as 



Foil) 



cosh(V2) 



sinh(//2) -|-acosh(?/2) 



Pd - 2kTt^a^j^Foil) + kTm^j^F^,aw)(.l), (63) 
where I = Lm, a — m^, a = with 7 = 7 = /i*/2/i. Typical values for these parameters are 

kT = 4. 10'^^ J, fJ.^ 0.2 mM ^ m - 0.05 nm"\ a = a 2.0 . 
The coefficients in Eq. (^3|) are then 

Co = 2kBT^io? = 4000 Pa, ci = keTm^ ^ 500 Pa, e/eq ^ 80 , 

where the ratio e/eo is for water to air. The disjoining pressure Pd is given in Pascals. For these values a plot of P^ 
is shown in Fig. ^ The solid curve is the total value of Pd while the dashed curve gives the repulsive contribution 
from Pg and the attractive contribution from F(^ydW) is shown as the dotted curve. The theory predicts a collapse 
certainly by L = 8 nm and since the collapse corresponds to a first order phase transition the Maxwell construction 
will predict that collapse will be observed at larger L. Since we do not have a theory for the short-range repulsive 
force that eventually stabilizes the NBF we cannot use the Maxwell construction to give an accurate value for where 
the film becomes metastable but 8 nm < L < 20 nm would be a reasonable range. 

The attractive van der Waals contribution pjj"'''^^ above can be compared with the Hamaker form at short distances 
predicted by Eq. (|6^). However, numerical study of the case of interest here shows that the values of L for which this 
behavior holds are too small to be relevant to the collapse described above. 
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For P 3> 1/m the leading term is given by Eq. ( 61 ) and the value of a is crucial in determining the overall coefficient. 



Indeed, for a ^ 1 the behavior will given by non-leading terms not shown here. In Fig. ^ we plot p^^"^^"^ versus L. It 
is seen that while the large L asymptotic form ( |6l| ) is a good approximation for L > 1/m = 20nm, in the the region 
important to the collapse, 5nm < L < 20nm, the full result deviates strongly from this form. Thus we find that the 
full expression for pj"''^^ must be used in the region of interest. 

The surface-charging mechanism is very important to the prediction of the collapse transition. If a fixed surface 
charge is used we should omit the quadratic term in <p in the expansion of the source in Eqs. (^^5|). This corresponds 
to setting a — where is occurs explicitly in the expressions for Pq and F^^dw) in Eq. ( ]6^ ) while not changing the 
values of the coefficients cq, ci. In this case for the parameters above there is no collapse. 

In a more general case with both /i!J^ , //I non-zero the same expressions as in Eq. (^3|) apply but with the gen- 
eralization a ^ a. The effect of choosing fixed a = 2.0 but varying a can be seen in Fig. p where the values used 
above of cq = 4000 Pa, c\ = 500 Pa are adopted but different values of a = 1.0, 1.5,2.0,3.0 are used. The term 
most affected is the repulsive mean-field term and even in the linearized theory this is very sensitive, as we should 
expect, to the charging mechanism for the surfaces. We see that the predicted properties of the collapse transition 
are strongly dependent on the choice of a and a and therefore on the details of the film being studied. 
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VII. BEHAVIOR OF THE SURFACE CHARGE 



In the limit 5-^0 the surface charge (on one surface) per unit area, cr, can be seen to be 



d 



d 



e / 5 d 



/3J, 



(64) 



where 7, 7 are defined in Eq. (p^. 

Even in the Gaussian approach the source terms encode the non-trivial charging properties of the surfaces bounding 
the ionic liquid. From Eq. (^) we have used the approximation for the source 



/i*^ exp(i/3e0) + exp(— i/Jet; 



iX(j) - }^l3einn?-i(f + 



(65) 



The term linear in (f) represents a fixed surface charge but the term in (jx^ corresponds to surface charge fluctuations. 

We decompose a in terms of the mean field contribution and the van der Waals contribution cr = ctq + Uvdw ■ The 
term coming from the mean field (and constant or ideal) contributions to the grand potential is 



ao(i) = 



m (a + tanh(2|i^)) 



, ,mL. 
2tanh( — ) 



(a + tanh(ifi)) 



The contribution from the van der Waals term is independent of 7 and can be written as 



We find 



2tt 



dk k 



CTvdW 



fik) 



^- ^ RJ 

2 da 



where / is given by 



1-P{k) (fceoA + a + yFTT)2 
fceo/e + a - \/I?M~T 



exp 



(-/v/pTT), 



keo/t 



exp 



{~l\/k^ 



(66) 



(67) 



(68) 



(69) 



For the parameters given in the previous section: n — 0.2 mM, m = 0.05 nm~^, a = a = 2.0 the behavior of a is 
shown versus L in Fig. ^. One notices that the surface charge is regularized on varying L, though not drastically 
(about 10% over 40 nm) just up to the film thickness L ~ lOnm. One sees that in this case the effect of the van der 
Waals term is to decrease the value of cr from its mean field value. 



VIII. CONCLUSIONS AND OUTLOOK 



In this paper we have presented a field theoretic formulation of the electrostatic interactions in soap film like 
systems, which treats on the same footing the zero frequency van der Waals or Lifshitz terms and the contributions 
coming from ionic fiuctuations. The basic idea is to use the static part of the QED Lagrangian coupled to the charge 
density coming from the ions in the system and then integrating over the electrostatic potential and the positions of 
the ions (which are treated classically) . The time dependent and magnetic field terms in the full QED Lagrangian are 
thus neglected, this is equivalent to the non retarded limit where the velocity of light c ^ 00. Retardation effects can 
be taken into account by summing over the non zero Matsubara frequencies, however the coupling of these terms with 
the ionic distribution is weak. The incorporation of retardation effects requires in addition the frequency dependence 
of the electric permittivies. 

This treatment is easily applicable to systems with spatially varying dielectric constants and elegantly avoids 
calculations of the arising image effects. The formalism also allows the incorporation of surface charges induced by 
equilibrium processes. In the grand canonical ensemble one obtains a Sine-Gordon field theory. Linearizing this 
theory leads to a soluble Gaussian field theory and is equivalent to the Debye-Hiickel approximation, which should 
be valid for weak ionic concentrations. In this form the evaluation of the grand potential is carried out by using the 
Feynman kernel for simple harmonic oscillators. The use of the Feynman kernel in the field theoretic formalism allows 
us reproduce a wide range of results established in the literature via other methods. The effect of surface charge 
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fluctuations are can also be taken into account and it was also shown how many different layers of varying dielectric 
constant simply lead to a renormalization of the simple bilayer result. 

Preliminary investigation of the relevant experimental soap film parameters shows that the van der Waals contri- 
butions lead to a weak (screened) attraction at large intersurface separations. For thinner films there is an increased 
attraction, which can overcome the mean field repulsion present in the models considered here. However in the region 
where the film collapses, the attraction does not have the simple Hamaker form and is strongly dependent on 

the Debye mass m and the surface surface-charging parameters. Indeed an essential ingredient is the inclusion of the 
surface-charging mechanism which, although treated here in the Gaussian approximation, nevertheless predicts an 
L-dependent surface charge density which is important to the details of the collapse. A linear approximation to the 
surface charge source is inadequate since it leads to a fixed surface charge and hence a diverging mean-field repulsive 
pressure Pq as L decreases. By including the term in (fp' in the expansion of the sources, the surface charge is shown 
to decrease as L decreases and so the divergence in the mean-field repulsion Pq is regulated. 

Although not necessary for the theory presented here it is interesting to compare the outcome with our work on 
the ID coulomb gas model for a soap film ||l9|. The mechanism in ID for the collapse was the changing balance of 
contributions to Pd between the even and odd eigenfunctions of the Mathieu equation as the film thickness L varied. 
The important states were the lowest lying ones including the ground-state. In the Gaussian approximation used 
in the present work we can ask which are the important eigenfunctionals of the theory in 3D which play a similar 
role. In this model the eigenfunctionals are products of harmonic oscillator eigenfunctions for each of the transverse 
momentum modes separately for which the coordinate is (j){p). The source term is, however, only a function of the 
zero mode ^(0) and the important corresponding term in the wavefunctional from Eq. ( p^ is (p = 0, X = (^(0)) 

= /i„(X)exp(-mXV2) , 

where n is the oscillator excitation number. For large n the Hermite polynomial has an oscillatory factor 

hn{X) ^ COS {\/n mX — {n — \)tt / 2) . 

The relevant coefficient is the overlap of this wavefunction with the source term in Eq. (|l]) exp(i\/]4/^*/3eX) and it 
is clear that this will not be large unless the oscillatory factors match. This will be the case when 

\/n m ^ vC4/i*/3e , 

and the corresponding energy values are of order E '-^ (n-\- l/2)m ~ Afji*'^{[3ef' /m giving an extensive contribution to 
the free energy as we must expect. Hence the important states are not the ground-state and those nearby but highly 
excited states which carry the extensive nature of the system. 

By taking experimentally reasonable values for the parameters in our formulae we obtain acceptable predictions 
for the collapse phenomenon and surface-charging for simple description of the film. However, to make accurate 
predictions will require the film to be modeled as a multiple layer with the correct permittivities for each layer 
and possible charging potentials included. We must also include the contributions of non-linear and non-Gaussian 
operators and to do this involves three ingredients. The first is to solve the non-linear mean-field equations, the 
second is to develop the perturbation theory for the non-Gaussian source operators within the Gaussian field theory 
and the third to use perturbation theory for the non-Gaussian interactions given by the Sine-Gordon theory in the 
film interior. Because the system is not translation invariant these perturbation theories is not standard but it has 
been developed and will be presented in a succeeding paper together with the non-linear mean-field formulae. 
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FIG. 1: The disjoining pressure Pd for the linearized theory using typical parameter values: /iNaCi = 0.2 mM, moebye = 
0.05 nm"^ , a = 2.0. The solid line is the full result, the dashed line is the repulsive contribution and the dotted line the 
attractive, van der Waals contribution. 
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FIG. 2: The van der Waals contribution P^" for the hnearized theory plotted against L using typical parameter values: 
MNaCi = 0.2 mM, moebye = 0.05 nm~^ , a = 2.0. The dashed line is the asymptotic formula Eq. (pm which is seen to be 
accurate until L < 1/m — 20nm. The standard Hamaker form oc L"^ is not applicable for the relevant values of L. 
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FIG. 3: The disjoining pressure Pd for the linearized theory using typical parameter values: /XNaCi ~ 0.2 mM, mDobyc ~ 
0.05 nm~^ , a — 2.0 but for different values of a occurring explicitly in Eq. (|63|). The curves from top down are for 
a = 0.5, 1.0, 1.5,2.0,3.0. 
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FIG. 4: The surface charge a in rnC plotted versus L for /LiNaCi = 0.2 mM, moebye = 0.05 nm~ , a = 2.0. The sohd hne is 
the full result, the dashed line is the mean field contribution uo. The van der Waals contribution, cr(„dvv) is negligible until 
L < lOnm where it is a responsible for the rapid decrease in a. 
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